function MergeNatureResults(n1)

ValueFolderTemp='/sps/isc/rneveu/EtudeIRM/Nature/ResultsJournal/';
VecA=192:1:198;
VecA=[VecA 65];
Veca=224:1:230;
Veca=[Veca 97];
VecC=[199 67];
Vecc=[231 99];
VecE=200:203;
VecE=[VecE 69];
Vece=232:1:235;
Vece=[Vece 101];
VecI=204:207;
VecI=[VecI 73];
Veci=236:1:239;
Veci=[Veci 105];
Vecn=[241 110];
VecO=210:214;
VecO=[VecO 216 79];
Veco=242:1:246;
Veco=[Veco 248 111];
VecU=217:1:220;
VecU=[VecU 85];
Vecu=249:1:252;
Vecu=[Vecu 117];
VecY=[221 89];
Vecy=[253 121];
Vecs=[223 115];
VecSpecialCharacters=[{VecA} {Veca} {VecC} {Vecc} {VecE} {Vece} {VecI} {Veci} {Vecn} {VecO} {Veco} {VecU} {Vecu} {VecY} {Vecy} {Vecs}];
VecValueDateTemp={'Jan' 'Feb' 'Mar' 'Apr' 'May' 'Jun' 'Jul' 'Aug' 'Sep' 'Oct' 'Nov' 'Dec'};
ValueClock=clock;
disp(['Start of the processing at ' num2str(ValueClock(1,4)) ':' num2str(ValueClock(1,5)) ':' num2str(ValueClock(1,6)) ' for journal ' num2str(n1)]);
ValueFile=['ResultsJournal' num2str(n1) '.mat'];
load([ValueFolderTemp ValueFile]);
n3=2020;
eval(sprintf('VecTemp=VecResultsArticle%d;',n3));
n3=n3-1;
VecResults=VecTemp;
try
    load([ValueFolderTemp 'ResultsJournalTemp2Temp' num2str(n1) '.mat']);
end
while n3>=1860
    n3
    eval(sprintf('VecTemp=VecResultsArticle%d;',n3));
    for n5=1:size(VecTemp,1)
        n4=1;
        TestTemp=0;
        ValueAuthor=VecTemp{n5,1};
        while n4<=size(VecResults,1)&&TestTemp==0
            ValueTemp2=VecResults{n4,1};
            if strcmps(ValueTemp2,ValueAuthor,VecSpecialCharacters)>=0.99
                ValueTemp4=VecResults{n4,2};
                ValueTemp5=VecTemp{n5,2};
                ValueTemp4=[ValueTemp4;ValueTemp5];
                VecResults(n4,2)={ValueTemp4};
                ValueTemp4=VecResults{n4,3};
                ValueTemp5=VecTemp{n5,3};
                ValueTemp4=[ValueTemp4;ValueTemp5];
                VecResults(n4,3)={ValueTemp4};
                ValueTemp4=VecResults{n4,4};
                ValueTemp5=VecTemp{n5,4};
                ValueTemp4=[ValueTemp4;ValueTemp5];
                VecResults(n4,4)={ValueTemp4};
                ValueTemp4=VecResults{n4,5};
                ValueTemp5=VecTemp{n5,5};
                ValueTemp4=[ValueTemp4;ValueTemp5];
                VecResults(n4,5)={ValueTemp4};
                ValueTemp4=VecResults{n4,6};
                ValueTemp5=VecTemp{n5,6};
                ValueTemp4=[ValueTemp4;ValueTemp5];
                VecResults(n4,6)={ValueTemp4};
                ValueTemp4=VecResults{n4,7};
                ValueTemp5=VecTemp{n5,7};
                ValueTemp4=[ValueTemp4;ValueTemp5];
                VecResults(n4,7)={ValueTemp4};
                TestTemp=1;
            end
            n4=n4+1;
        end
        if TestTemp==0
            ValueTemp6=size(VecResults,1)+1;
            VecResults(ValueTemp6,1)={ValueAuthor};
            ValueTemp4=VecTemp{n5,2};
            VecResults(ValueTemp6,2)={ValueTemp4};
            ValueTemp4=VecTemp{n5,3};
            VecResults(ValueTemp6,3)={ValueTemp4};
            ValueTemp4=VecTemp{n5,4};
            VecResults(ValueTemp6,4)={ValueTemp4};
            ValueTemp4=VecTemp{n5,5};
            VecResults(ValueTemp6,5)={ValueTemp4};
            ValueTemp4=VecTemp{n5,6};
            VecResults(ValueTemp6,6)={ValueTemp4};
            ValueTemp4=VecTemp{n5,7};
            VecResults(ValueTemp6,7)={ValueTemp4};
        end
    end
    n3=n3-1;
    if mod(n3,50)==0
        save([ValueFolderTemp 'ResultsJournalTemp2Temp' num2str(n1) '.mat']);
    end
end
for n3=1:size(VecResults,1)
    ValueTemp=VecResults{n3,7};
    ValueResultsTemp2=VecResults{n3,2};
%         if size(ValueResultsTemp2,1)==0
%             ValueResultsTemp2=[NaN NaN];
%         end
    ValueResultsTemp3=VecResults{n3,3};
%         if size(ValueResultsTemp3,1)==0
%             ValueResultsTemp3=[1 1 0 0 0 1];
%         end
    ValueResultsTemp4=VecResults{n3,4};
%         if size(ValueResultsTemp4,1)==0
%             ValueResultsTemp4=zeros(1,6);
%         end
    ValueResultsTemp5=VecResults{n3,5};
%         if size(ValueResultsTemp5,1)==0
%             ValueResultsTemp5=zeros(1,6);
%         end
    ValueResultsTemp6=VecResults{n3,6};
%         if size(ValueResultsTemp6,1)==0
%             ValueResultsTemp6={''};
%         end
    if size(ValueTemp,1)>=1&&size(ValueResultsTemp2,1)>=1&&size(ValueResultsTemp3,1)>=1&&size(ValueResultsTemp4,1)>=1&&size(ValueResultsTemp5,1)>=1&&size(ValueResultsTemp6,1)>=1
        VecDate=[];
        for n4=1:size(ValueTemp,1)
            ValueDate=ValueTemp{n4,1};
            ValueYear=ValueDate{1,1};
            ValueMonth=ValueDate{1,2};
            ValueDay=ValueDate{1,3};
            if size(ValueYear,2)>=1
                ValueYearTemp=str2num(ValueYear);
            else
                ValueYearTemp=NaN;
            end
            if size(ValueMonth,2)>=1
                ValueMonthTemp=str2num(ValueMonth);
                if size(ValueMonthTemp,2)==0
                    n61=1;
                    TestTemp3=0;
                    while n61<=12&&TestTemp3==0
                        if strcmpi(ValueMonth,VecValueDateTemp{1,n61})==1
                            ValueMonthTemp=n61;
                            TestTemp3=1;
                        end
                        n61=n61+1;
                    end
                end
            else
                ValueMonthTemp=NaN;
            end
            if size(ValueDay,2)>=1
                ValueDayTemp=str2num(ValueDay);
            else
                ValueDayTemp=NaN;
            end
            VecDate=[VecDate;ValueYearTemp ValueMonthTemp ValueDayTemp];
        end
        [ValueTemp2,ValueTemp3]=sort(VecDate(:,1));
        VecDateTemp=[ValueTemp2 VecDate(ValueTemp3,2) VecDate(ValueTemp3,3)];
        ValueResultsTemp2=ValueResultsTemp2(ValueTemp3,:);
        ValueResultsTemp3=ValueResultsTemp3(ValueTemp3,:);
        ValueResultsTemp4=ValueResultsTemp4(ValueTemp3,:);
        ValueResultsTemp5=ValueResultsTemp5(ValueTemp3,:);
        ValueResultsTemp6=ValueResultsTemp6(ValueTemp3,:);
        n5=1;
        for n4=1:size(VecDateTemp,1)-1
            if VecDateTemp(n4+1,1)>VecDateTemp(n4,1)&&n4>n5
                n6=n4;
                ValueTemp2=VecDateTemp(n5:n6,:);
                [ValueTemp3,ValueTemp4]=sort(ValueTemp2(:,2));
                ValueTemp2=[ValueTemp2(ValueTemp4,1) ValueTemp3 ValueTemp2(ValueTemp4,3)];
                VecDateTemp(n5:n6,:)=ValueTemp2;
                ValueResultsTemp22=ValueResultsTemp2(n5:n6,:);
                ValueResultsTemp22=ValueResultsTemp22(ValueTemp4,:);
                ValueResultsTemp2(n5:n6,:)=ValueResultsTemp22;
                ValueResultsTemp22=ValueResultsTemp3(n5:n6,:);
                ValueResultsTemp22=ValueResultsTemp22(ValueTemp4,:);
                ValueResultsTemp3(n5:n6,:)=ValueResultsTemp22;
                ValueResultsTemp22=ValueResultsTemp4(n5:n6,:);
                ValueResultsTemp22=ValueResultsTemp22(ValueTemp4,:);
                ValueResultsTemp4(n5:n6,:)=ValueResultsTemp22;
                ValueResultsTemp22=ValueResultsTemp5(n5:n6,:);
                ValueResultsTemp22=ValueResultsTemp22(ValueTemp4,:);
                ValueResultsTemp5(n5:n6,:)=ValueResultsTemp22;
                ValueResultsTemp22=ValueResultsTemp6(n5:n6,:);
                ValueResultsTemp22=ValueResultsTemp22(ValueTemp4,:);
                ValueResultsTemp6(n5:n6,:)=ValueResultsTemp22;
                n5=n4+1;
            end
        end
        if n4>n5
            n6=n4;
            ValueTemp2=VecDateTemp(n5:n6,:);
            [ValueTemp3,ValueTemp4]=sort(ValueTemp2(:,2));
            ValueTemp2=[ValueTemp2(ValueTemp4,1) ValueTemp3 ValueTemp2(ValueTemp4,3)];
            VecDateTemp(n5:n6,:)=ValueTemp2;
            ValueResultsTemp22=ValueResultsTemp2(n5:n6,:);
            ValueResultsTemp22=ValueResultsTemp22(ValueTemp4,:);
            ValueResultsTemp2(n5:n6,:)=ValueResultsTemp22;
            ValueResultsTemp22=ValueResultsTemp3(n5:n6,:);
            ValueResultsTemp22=ValueResultsTemp22(ValueTemp4,:);
            ValueResultsTemp3(n5:n6,:)=ValueResultsTemp22;
            ValueResultsTemp22=ValueResultsTemp4(n5:n6,:);
            ValueResultsTemp22=ValueResultsTemp22(ValueTemp4,:);
            ValueResultsTemp4(n5:n6,:)=ValueResultsTemp22;
            ValueResultsTemp22=ValueResultsTemp5(n5:n6,:);
            ValueResultsTemp22=ValueResultsTemp22(ValueTemp4,:);
            ValueResultsTemp5(n5:n6,:)=ValueResultsTemp22;
            ValueResultsTemp22=ValueResultsTemp6(n5:n6,:);
            ValueResultsTemp22=ValueResultsTemp22(ValueTemp4,:);
            ValueResultsTemp6(n5:n6,:)=ValueResultsTemp22;
        end
        n5=1;
        for n4=1:size(VecDateTemp,1)-1
            if VecDateTemp(n4+1,2)~=VecDateTemp(n4,1)&&n4>n5 % Error here but without any consequence over the following steps as the find function is used for both years and months. The original goal of sorting all articles by increasing order of date was useless
                n6=n4;
                ValueTemp2=VecDateTemp(n5:n6,:);
                [ValueTemp3,ValueTemp4]=sort(ValueTemp2(:,3));
                ValueTemp2=[ValueTemp2(ValueTemp4,1) ValueTemp2(ValueTemp4,2) ValueTemp3];
                VecDateTemp(n5:n6,:)=ValueTemp2;
                ValueResultsTemp22=ValueResultsTemp2(n5:n6,:);
                ValueResultsTemp22=ValueResultsTemp22(ValueTemp4,:);
                ValueResultsTemp2(n5:n6,:)=ValueResultsTemp22;
                ValueResultsTemp22=ValueResultsTemp3(n5:n6,:);
                ValueResultsTemp22=ValueResultsTemp22(ValueTemp4,:);
                ValueResultsTemp3(n5:n6,:)=ValueResultsTemp22;
                ValueResultsTemp22=ValueResultsTemp4(n5:n6,:);
                ValueResultsTemp22=ValueResultsTemp22(ValueTemp4,:);
                ValueResultsTemp4(n5:n6,:)=ValueResultsTemp22;
                ValueResultsTemp22=ValueResultsTemp5(n5:n6,:);
                ValueResultsTemp22=ValueResultsTemp22(ValueTemp4,:);
                ValueResultsTemp5(n5:n6,:)=ValueResultsTemp22;
                ValueResultsTemp22=ValueResultsTemp6(n5:n6,:);
                ValueResultsTemp22=ValueResultsTemp22(ValueTemp4,:);
                ValueResultsTemp6(n5:n6,:)=ValueResultsTemp22;
                n5=n4+1;
            end
        end
        if n4>n5
            n6=n4;
            ValueTemp2=VecDateTemp(n5:n6,:);
            [ValueTemp3,ValueTemp4]=sort(ValueTemp2(:,3));
            ValueTemp2=[ValueTemp2(ValueTemp4,1) ValueTemp2(ValueTemp4,2) ValueTemp3];
            VecDateTemp(n5:n6,:)=ValueTemp2;
            ValueResultsTemp22=ValueResultsTemp2(n5:n6,:);
            ValueResultsTemp22=ValueResultsTemp22(ValueTemp4,:);
            ValueResultsTemp2(n5:n6,:)=ValueResultsTemp22;
            ValueResultsTemp22=ValueResultsTemp3(n5:n6,:);
            ValueResultsTemp22=ValueResultsTemp22(ValueTemp4,:);
            ValueResultsTemp3(n5:n6,:)=ValueResultsTemp22;
            ValueResultsTemp22=ValueResultsTemp4(n5:n6,:);
            ValueResultsTemp22=ValueResultsTemp22(ValueTemp4,:);
            ValueResultsTemp4(n5:n6,:)=ValueResultsTemp22;
            ValueResultsTemp22=ValueResultsTemp5(n5:n6,:);
            ValueResultsTemp22=ValueResultsTemp22(ValueTemp4,:);
            ValueResultsTemp5(n5:n6,:)=ValueResultsTemp22;
            ValueResultsTemp22=ValueResultsTemp6(n5:n6,:);
            ValueResultsTemp22=ValueResultsTemp22(ValueTemp4,:);
            ValueResultsTemp6(n5:n6,:)=ValueResultsTemp22;
        end
        VecResults(n3,2)={ValueResultsTemp2};
        VecResults(n3,3)={ValueResultsTemp3};
        VecResults(n3,4)={ValueResultsTemp4};
        VecResults(n3,5)={ValueResultsTemp5};
        VecResults(n3,6)={ValueResultsTemp6};
        VecResults(n3,7)={VecDateTemp};
    end
end
save([ValueFolderTemp 'ResultsJournalTemp' num2str(n1) '.mat'],'VecResults');
ValueClock=clock;
disp(['End of the processing at ' num2str(ValueClock(1,4)) ':' num2str(ValueClock(1,5)) ':' num2str(ValueClock(1,6)) ' for journal ' num2str(n1)]);